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Abstract:  Because  of  the  interest  in  predicting  the  characteristics 

of  laminar,  heated,  water  boundary  layers,  an  existing  code 
was  adapted  for  use  on  The  Pennsylvania  State  University's 
IBM  370  computer  and  extended  to  include  the  axisymmetric 
case.  The  similar  type  solutions  of  the  boundary  layer 
equations  made  it  possible  to  perform  a parametric 
study  and  to  use  the  results  in  developing  a means  for 
specifying  a stabilizing  temperature  distribution  on  a 
body.  Additional  studies  into  the  effect  of  temperature  on 
laminar  separation  and  transition  have  resulted  in  simplified 
criteria  for  predicting  these  phenomena. 
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Introduction 

Effective  calculation  of  laminar  boundary  layer  properties  has  been 
possible  either  by  the  approximate  Karman-Pohlhausen  method  or  by  more 
exact  numerical  Integration  methods.  See,  for  example,  References 
1,  2 and  3.  In  recent  years,  this  kind  of  calculation  has  been  extended 
to  the  case  of  the  heated,  water  boundary  layer.  In  Reference  4,  as 
part  of  the  computations  of  the  stability  of  heated,  laminar  boundary 
layers  In  water,  the  solution  of  the  mean-flow  equations  was  coded 
and  specific  results  were  presented.  Because  of  the  Interest  In  applying 
these  predictions  to  design  techniques,  this  part  of  the  code  was  adapted 
to  The  Pennsylvania  State  University's  IBM  370  computer.  The  code, 
as  It  appeared  In  Reference  4,  treated  the  two-dimensional  flow  case 
and  therefore  had  to  be  altered  for  axisymmetrlc  flow.  This  memorandum 
summarizes  the  development  of  the  altered  code  and  the  results  of  calcula- 
tions using  the  code  to  date.  It  also  summarizes  some  of  the  uses 
to  which  the  output  data  from  the  code  have  been  put.  This  includes 
the  development  of  a procedure  to  predict  the  temperature  distribution 
necessary  to  insure  stability  of  the  laminar  boundary  layer. 


Development  of  Equations 


The  basic  boundary  layer  equations  were  derived  for  an  incompressible. 


axisymmetrlc  flow  with  heat  addition  at  the  boundary.  Such  effects  as 
buoyancy  and  dissipation  of  energy  by  friction  were  ignored.  The  resulting 
equations,  including  an  equation  of  energy  balance  are  as  follows: 
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9(pru)  . 3 (pry)  _ „ 

3x  3y  ^ 


, 3u  3u.  *^^e  3 , Su. 

^ ^ 3^  ay-> 


^ . 3T  . 3T.  ^ .3  3T, 

P 3^  ^ 37^  = ^ ^ 37^ 


There  is  another  momentum  equation  in  the  z direction  which  is  not  coupled 
with  the  above  and,  of  course,  in  the  y direction,  3p/3y  = 0.  The  x and  y 
coordinates  are  curvilinear  coordinates  tangent  to  and  perpendicular  to 
the  body  surface  respectively.  The  axisjnnmetrlc  nature  of  the  flow  is 
reflected  in  the  continuity  equation  with  the  inclusion  of  the  r term. 

This  represents  a departure  from  Lowell  and  Reshotko  in  Reference  4. 

In  order  to  solve  these  equations,  they  are  transformed  to  ordinary 
differential  equations  by  defining  a stream  function,  if),  and  a nondimens ional 
y coordinate,  n. 


'I'  = 


n = (u  /v 

' ' e ~ 


p dy 


(2) 


Making  use  of  these  definitions  and  Equations  (1),  the  following  final 
equations  result: 


• j.  M/—  _ j.  /M  + X + 1^  _ Q 


(p  p £")•  + M(|  - f^)  + (•  2 


(3) 


PC  (^-  " f0'  + (pk0')' 

P 2 
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The  primes  indicate  derivatives  with  respect  to  n*  These  equations  are 
explicitly  independent  of  x and  the  velocity  and  temperature  profiles 
that  result  are  given  by  similarity  solutions.  0 is  a normalized 
temperature  difference  and  M and  X are  given  by: 


and 


M 


3c 

U 


= 2x  ^ 
r dx 


(4) 


The  following  boundary  conditions  apply  for  Equations  (1) : 


aty  = 0 u = v=  0 

T = T 

w 


as  y 00 


u U 

e 

00 


(5) 


In  terms  of  the  transformed  variables  in  Equations  (3),  these  become: 


atn=0  f'=f=0 

0 = 1 

as  n " f'  1 

0 -►  0 


(6) 


The  solution  process  requires  that  the  physical  properties  of  the  water 
(Cp,  p,  M,  k)  be  known  at  each  location  in  the  boundary  layer.  Since  the 
temperature  varies  through  the  boundary  layer,  so  will  the  physical 
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properties  and  one  must  know  their  values  and  the  values  of  their 
derivatives  with  temperature.  The  formulas  (from  Reference  4)  for 
the  variation  of  the  water  properties  with  temperature  are  presented 
in  Appendix  A.  The  derivatives  are  obtained  directly  from  these  equations. 

Equations  (3)  are  nonlinear  and  are  coupled  through  the  fluid 
properties.  In  addition,  there  are  problems  of  initial  conditions 
(f"  and  0'  at  the  wall)  and  of  the  value  of  r)  at  which  to  terminate  the 
integration.  The  method  of  Nachtscheim  and  Swigert  (Reference  5)  was 
adopted  for  the  numerical  integration  process  and  was  found  to  be 
satisfactory.  All  this  is  discussed  in  detail  in  Reference  4. 

Boundary  Layer  Computer  Program 

The  computer  program  for  the  boundary  layer  calculations  was  extracted 
from  the  more  general  code  of  Reference  4 and  adapted  to  the  University's 
IBM  370  computer.  Besides  being  altered  to  account  for  computer 
Idiosyncracies,  changes  were  made  to  Include  the  axisymmetric  body  case 
and  to  calculate  additional  boundary  layer  quantities  that  were  not 
available  in  the  original  program.  The  inclusion  of  the  axisymmetric 
case  results  in  the  term  X appearing  in  Equations  (3).  For  the  two- 
dimensional  case,  X is  zero. 

The  total  code  consists  of  a Main  program  and  a Mean-Flow  Equation 
(MFEQN)  program  plus  three  additional  subroutines.  The  Main  and  MFEQN 
programs  could  have  been  combined,  but  were  left  separate  in  deference 
to  the  original  makeup  of  the  code.  The  code  originally  included 
stability  calculations  and  the  Main  program  was  used  to  control  overall 
execution  and  as  such,  was  used  in  conjunction  with  eleven  subroutines. 

The  MFEQN  routine  controls  the  solution  of  the  mean-flow  equations  and 
to  do  this  calls  upon  the  DIFF,  TEMVAR,  and  ADAMS  subroutines.  The 
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DIFF  routine  evaluates  the  first  order  differential  mean-flow  equations 
and  calls  upon  TEMVAR  to  evaluate  required  mean-flow  properties  and 
their  derivatives.  The  mean-flow  equations  are  integrated  in  double 
precision  by  the  ADAMS  routine.  The  integration  routine  is  designed  to 
use  single  step  Runge-Kutta  procedures  alone  or  in  combination  with  the 
multi-step  Adams-Moulton  predictor-corrector  scheme. 

The  key  quantities  that  are  needed  as  Inputs  to  the  boundary  layer 
program  are  the  free-stream  temperature,  the  wall  temperature,  the 
Falkner-Skan  parameter,  M,  and  the  radius  gradient  parameter,  X.  The 
M and  X values  could  have  been  input  directly  or  been  calculated  from  the 
provided  quantities  that  make  up  these  parameters  (l.e.,  x,  U^,  dU^/dx, 
r,  dr/dx).  The  latter  quantities  were  the  ones  that  were  used  in 
practice.  By  their  use,  actual  values  of  the  various  boundary  layer 
thicknesses,  Reynolds  numbers,  etc.,  could  be  calculated. 

The  important  output  quantities  Include  a term  referred  to  as 
ri^*,  the  shape  factor,  H=6  /O,  the  wall  shear  stress,  T,  and  the  heat 

flux,  q.  The  term  Hg*  is  derived  in  the  evaluation  of  the  displacement 

* 

thickness  o . It  should  be  mentioned  at  this  time  that  both  the 
displacement  and  momentum  thicknesses  are  defined  in  the  sense  of  two- 
dimensional  flow  definitions.  This  leads  eventually  to  slight  differences 
when  comparing  with  the  same  quantities  defined  differently  in  other 
sources.  Little  difference  in  the  shape  factor  should  result,  however, 
when  the  point  on  the  body  under  consideration  has  a small  boundary 
layer  thickness  compared  to  the  body  radius. 

The  displacement  thickness  for  the  heated  water  boundary  layer 
is  here  used  to  represent  the  decrease  in  mass  flow: 
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P»  U 6 

00  g 


w 

1 "e 


-pu)dy 


(1  - p u)dy 


This  can  be  transformed  in  terms  of  dn,  so  that 


f (i  - u, 

• o'’ 


The  integral  portion  of  this  equation  is  what  is  referred  to  as  Hr*- 

o 

Therefore, 


5 ^ U ^ *^6* 


In  a similar  manner,  the  momentum  thickness  is  defined  as 


I p u 


(l-u)dy  , 


or  in  terms  of  dri 


e = ( y 

e 


The  shape  factor  is,  therefore: 


j U(l-u) 


H ••  TT-  “ 

6 “ 


I u(l-u) 


(13) 
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Results  of  Boundary  Layer  Computations 

The  code  was  used  primarily  to  run  a sufficient  number  of  cases 
of  varied  parameters  so  that  one  could  generate  design  information 
and  procedures  for  heated  axisymmetric  bodies.  Checks  were  first  made 
with  the  results  presented  in  Reference  4.  Much  of  the  information 
shown  in  Reference  4 is  in  the  form  of  plots  which  are  difficult  to 
read  accurately.  Some  of  the  quantities  which  are  listed  in  tables  in 
Reference  4 and  which  can  be  compared  are  those  listed  in  Table  1. 

It  is  apparent  that  the  comparisons  are  quite  good  and  the  largest 
difference  is  of  the  order  of  2 in  the  sixth  significant  figure.  Since 
these  are  key  terms  in  the  scheme  of  calculations,  it  is  assumed  that 
all  the  other  calculated  quantities  are  of  a similar  comparative  accuracy. 

Another  comparison  (Table  2)  is  made  between  the  results  of  the 
unheatec  case  with  no  pressure  gradient  and  the  classical  Blasius  solution. 
Although  the  displacement  thickness  results  appear  to  check  very  closely, 
the  momentum  thickness  results  appear  to  differ  by  about  0.5%.  The 
results  were  considered  close  enough  to  verify  the  operation  of  the 
present  computer  program. 

The  bulk  of  the  boundary  layer  computations  were  performed  as  a 
parametric  study  in  which  M,  X,  and  AT  were  the  parameters.  In  all 
cases  the  free-stream  temperature,  T^,  was  kept  at  bCF.  Also,  so  that 
comparative  values  of  shear  stress  and  heat  flux  could  be  obtained, 
a value  of  arc  length,  x,  was  made  equal  to  4 ft  and  the  velocity  at  the 
edge  of  the  boundary  layer,  U^,  was  chosen  to  be  50  ft/sec.  A summary 
of  the  results  of  this  parameter  study  is  to  be  found  in  Table  3.  Other 
computed  quantities  could  have  been  included  in  this  table,  but  it  was 
thought  that  these  are  the  more  important  ones  for  later  use.  Plots  of 
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these  data  are  shown  in  Figures  1-15.  As  can  be  seen,  within  certain 
ranges  of  parameters,  the  data  are  fairly  abundant  and  amenable  to  inter- 
polation. The  ranges  in  which  seemingly  sparse  data  exist  were  generally 
found  to  show  up  Infrequently  in  practical  cases  or  were  beyond  the 
computer  program's  capabilities. 

An  additional  side  study  was  performed  in  order  to  learn  some  of  the 

general  effects  of  having  elevated  ambient  temperatures.  This  is  the  kind 

of  situation  that  would  occur  in  a water  tunnel  test  where  an  attempt  is 

being  made  to  extend  the  Reynolds  number  range  by  heating  the  tunnel 

water.  This  study  was  performed  for  the  two-dimensional  case  (X=0)  only. 

The  implication  here  is  that  at  least  a qualitative  idea  of  the  effect 

could  be  learned  and  serve  as  a guide  in  designing  an  experiment. 

Data  generated  for  a number  of  free-stream  temperatures  is  presented 

in  Table  4.  Plots  of  H,  q.  and  x are  shown  in  Figures  16-27.  Data 

for  a 60®F  free-stream  temperature  are  included  again  in  this  table  and 

these  figures  for  the  sake  of  comparison. 

Critical  Reynolds  Number  Correlation 

Stability  information  in  terms  of  the  critical  Reynolds  number,  R^*  , 

crit 

is  available  for  two-dimensional  flow  both  for  the  case  of  the  unheated 

boundary  with  pressure  gradient  and  for  the  flow  over  a flat  plate  with 

heating  at  the  wall.  The  results  of  the  unheated  case  are  found  in  Reference 

6 ‘^^e 

6 as  the  critical  Reynolds  number  versus  the  Pohlhausen  parameter  A = — , 

and  the  heated  case  results  were  obtained  from  Reference  4 as  critical 

Reynolds  number  versus  temperature  difference.  The  parameter  A and  AT 

* 

were  translated  into  the  shape  parameter  H = 6 /0  so  that  in  each 

case  Rg*  could  be  plotted  against  H.  The  curves  in  Figure  28  are 
°crlt 

the  result.  It  can  be  seen  that  the  correlation  is  quite  good,  indicating 
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that  the  stability  of  the  laminar  boundary  layer  is  strongly  dependent 

upon  H,  regardless  of  whether  it  is  obtained  by  favorable  pressure 

gradient  or  by  heat.  A similar  correlation  is  reported  in  Reference  7. 

This  correlation  provides  the  basis  for  the  development  of  the  type  of 

design  information  which  will  be  presented  here.  Although  based  on  two- 

dimensional  stability  information,  the  correlation  is  considered  valid 

for  the  axisymmetric  case  because  the  stability  equations,  under  the 

assumption  that  the  boundary  layer  thickness  is  small  compared  to  the  local 

body  radius,  are  the  same  for  both  cases. 

The  curve  in  Figure  28  for  the  variation  in  pressure  gradient  is 

repeated  in  Figure  29  and  extrapolated  in  the  low  R.*  range.  This 

*^crit 

curve  is  used  subsequently  to  represent  the  variation  of  H with  R^* 
regardless  of  how  the  H is  obtained.  Table  5 gives  specific  values 
from  this  curve  for  ease  in  effecting  computations  involving  this  curve. 
Determination  of  Body  Temperature  Distribution 

In  order  to  make  preliminary  design  estimates  of  the  temperatures 
needed  to  maintain  stable  flow  at  points  on  a body,  the  data  shown  in 
Table  3 were  used  In  conjunction  with  certain  criteria.  Essentially 
two  criteria  were  chosen  to  insure  the  maintenance  of  laminar  flow: 

(1)  the  provision  of  just  enough  heat  to  keep  the  Reynolds  number  (based 
on  displacement  thickness)  equal  to  the  critical  Reynolds  number  and, 

(2)  the  provision  of  enough  heat  so  that  the  peak  critical  Reynolds  number 
is  maintained.  These  are  referred  to  as  "minimum  heat"  and  "maximum 
heat"  conditions  respectively.  The  minimum  heat  criterion  implies 

that,  for  a particular  free-stream  velocity,  enough  heat  is  added  to 
make  the  operating  Reynolds  number,  R^*»  equal  to  the  critical  value. 

This  should  insure  that  there  will  never  be  any  amplification  of 
waves  in  the  laminar  boundary  layer.  The  maximum  heat  criterion  fixes 
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the  Rg*  at  its  maximum  value.  Whether  there  is  amplification  depends 
crit 

upon  the  free-stream  velocity  being  high  enough  to  have  the  operating 

R,*  exceed  the  maximum  R-* 

crit 

In  implementing  the  temperature  hunting  procedure,  the  data  in 

Table  3 can  be  filed  and  then  recovered  by  the  computer  for  interpolation 

purposes.  The  geometry  and  potential  flow  pressure  distribution  for 

a body  will  be  known  so  that  an  M and  a X can  be  determined  for  each 

point  on  the  body  under  consideration.  These  will  be  designated  as 

and  X^  for  a particular  body  point.  The  AT  required  at  a body  station 

can  be  determined  by  applying  one  of  the  criteria  already  mentioned 

and  interpolating  the  data  to  get  appropriate  quantities  corresponding 

to  M and  X . 
o o 

In  the  minimum  heat  case  the  procedure  would  be  as  follows: 

(1)  Determine  an  versus  AT  curve  corresponding  to  and  X^ 
as  illustrated  by  the  interpolation  procedure  pictured  in  Figure  30. 

The  curves  can  be  represented  by  spline  fits  and  values  are  extracted 
accordingly.  All  the  other  quantities  (n^*,  q,  t)  can  be  similarly 
handled  so  that  a curve  of  each  of  these  is  known  as  a function  of  AT 
for  particular  M^'s  and  X^'s. 

(2)  First  choose  AT=0.  The  corresponding  then  to  AT=0, 

A 

and  X^  can  be  used  to  determine  6 for  a particular  free-stream  velocity 
using  Equation  (10). 

(3)  Calculate  R^*  by  means  of 


h* 


* 


6 


U 

e 


This  is  the  operating  R^*. 


1 
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(4)  Enter  Figure  29  with  R^*  from  step  (3)  to  determine  a required 

H.  By  entering  with  the  operating  R^*  we  are  saying  that,  in  order  for 

this  to  be  the  R^*  , we  must  produce  a corresponding  value  of  H. 

crit 

(5)  The  H required  from  step  (4)  can  be  used  with  the  versus 

AT  curve  of  step  (1)  to  determine  the  required  AT.  This  required  AT 
now  will  change  the  value  of  originally  determined  in  step  (2)  for 

AT=0. 

(6)  Repeat  steps  (2),  (3),  (4),  and  (5)  until  the  AT  required 
converges  within  a desired  accuracy. 

(7)  Enter  the  q and  x versus  AT  curves  established  for  M and  X 

o o 

and  get  proper  q and  T values. 

(8)  Correct  the  q and  x values  for  proper  x and  values  according 
to  the  relationships: 


q = q(232.843)  /U  /x 

^corr.  ^ e 

X = x(0. 0056569)  U /U  /x 

corr.  e e 


(14) 


For  the  maximum  heat  case  the  somewhat  different  procedure  is  as 
follows: 

(1)  Assume  H=2.29.  This  corresponds  approximately  to  the  point 

where  R^*  reaches  its  maximum.  Adding  heat  beyond  the  value  that 
crit 

produces  the  maximum  R^*  will  result  in  R^*  becoming  smaller  (see 

°crlt  crit 

Reference  4)  and  thus  be  counterproductive. 

(2)  Enter  the  curve  of  versus  AT,  determined  in  step  (1)  of  the 
minimum  heat  procedure,  and  extract  the  AT  required  for  H^=2.29. 

(3)  Extract  q and  x for  the  AT  of  step  (2)  from  curves  of  q and 
X versus  AT  also  determined  in  step  (1)  of  the  minimum  heat  procedure. 
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Again,  the  same  corrections  for  and  x shovm  in  Equations  (14)  must 
be  made. 

In  addition  to  determining  the  local  values  of  q and  T,  one  can 
multiply  these  by  the  corresponding  local  area  elements  and  numerically 
sum  the  results  to  give  a running  total  of  heating  power  and  the  laminar 
skin  friction  drag  in  kilowatts  and  pounds  respectively.  The  formulas 
used  in  this  summation  are: 


N 

E 

n=0 


■ir(r  +r  . , ) (x  . ,-x  ) 
n n-fl^  n+1  n 

144 


2 X 10 


lam 


n=0 


144 


2 

n+1 


(15) 


Laminar  Separation 

An  attempt  was  made  to  obtain  some  limits  with  regard  to  the  M,  X, 
and  AT  parameters  beyond  which  laminar  separation  would  occur.  Using 
the  criterion  that  skin  friction  goes  to  zero  at  the  point  of  separation, 
curves  of  the  type  shown  in  Figures  8-11  were  extrapolated  to  obtain  the 
desired  limits.  Figure  31  shows  one  representation  of  limit  lines  so 
obtained.  If  the  operating  values  of  M and  X fall  below  and  to  the  left 
of  the  appropriate  AT  line,  laminar  separation  should  occur.  As  can 
be  seen,  temperature  difference  has  relatively  little  effect  on  laminar 
separation. 


Another  presentation  of  this  same  information  is  given  in  Figure  32. 
This  is  the  kind  of  plot  presented  in  Reference  8 for  the  two-dimensional 
(X=0.0)  flow  case.  The  curve  from  Reference  8 is  designated  by  "RAND" 
and  is  compared  with  the  X=0  case  from  the  present  study.  Although 
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showing  slightly  more  effect  of  temperature,  the  Rand  calculations 
also  indicate  the  minor  roll  of  heating  in  the  prevention  of  laminar 
separation. 

Transition 

The  prediction  of  transition  is  dependent  on  many  factors,  making 

it  difficult  to  develop  any  simplified  approach  to  such  predictions. 

Free-stream  turbulence  level,  surface  roughness,  body  vibrations,  etc., 

can  all  influence  the  location  of  transition.  In  addition,  body  shape, 

which  creates  the  body  pressure  distribution,  and  the  distribution  of 

heat  added  to  the  body  wall  can  materially  affect  the  location  of 

transition.  The  work  reported  here  concerns  only  the  effect  of  body 

shape  and  heat  addition  on  transition  and  ignores  the  other  effects. 

In  the  present  scheme  the  transition  location  can  be  found  without 

resorting  to  boundary  layer  or  stability  calculations.  Assuming  the 

body  shape  is  known,  one  would,  of  course,  have  to  generate  the  Falkner- 

Skan  and  radius-gradient  parameters. 

Through  the  courtesy  of  A.M.O.  Smith,  a plot  of  the  results  of 
9 

e stability  calculations  that  were  performed  on  a variety  of  heated 

and  unheated  wedges  was  obtained.  The  plot  appears  as  a band  of  data  when 

Rxtrans  plotted  against  values  of  shape  factor,  H.  (See 

U X 

Figure  33).  Rx..  is  defined  as  Rx  = — It  would  be 

logical  to  choose  the  lower  bound  of  this  band  as  the  transition  criterion. 
This  curve  is  reproduced  in  Figure  34  as  a plot  of  f^Xtj-ans^  versus 

H and  values  corresponding  to  points  on  the  curve  are  presented  in 
Table  6. 

Because  the  shape  factor  is  known  as  a function  of  the  parameters. 


M,  X,  and  AT  (see  Figures  1-4  and  Table  3),  the  transition  Reynolds 
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number  can  likewise  be  plotted  as  a function  of  these  parameters  as  shown 
in  Figures  35-38.  One  can  thus  interpolate  between  the  values  of  the 
parameters  and  know  what  the  transition  Reynolds  number  should  be. 

The  interpolation  can  be  performed  by  computer  using  procedures 
previously  outlined  and  the  curve  of  Figure  34. 

Some  transition  test  data  for  unheated  bodies  was  examined  and  a 
correlation  with  the  curve  of  Figure  34  was  attempted.  Most  of  these 
data  are  from  Reference  12  and  appear  in  Figure  34.  The  H values  for 
these  points  were  determined  from  the  data  of  Table  3.  Compared  with 
these  data,  it  is  obvious  that  the  curve  is  optimistic,  that  is,  it 
predicts  transition  at  a higher  Reynolds  number  and  thus,  at  a point 
farther  back  on  a body. 

The  prediction  of  the  AT  distribution  and  transition  location  for 
a specific  body  is  not  being  reported  in  this  memorandum.  This  will  be 
left  to  a later  memorandum  in  which  that  kind  of  result  will  be  correlated 
with  the  output  of  TAPS  (References  13-14) . Without  going  into  details 
of  this  correlation,  it  is  sufficient  to  say  that  the  H's  calculated 
by  TAPS  for  a given  temperature  distribution  have  been  consistantly 
higher  by  a factor  of  about  1.03,  than  those  used  in  obtaining  the 
temperature  distribution.  If  one  then  reduces  the  H's  employed  in  the 
determination  of  the  temperature  distribution  by  this  factor,  one  would 
expect  that  the  TAPS  calculations  with  the  new  temperatures  would 
produce  the  desired  H values.  Using  this  idea  in  reverse,  if  the  H 
values  of  the  data  points  in  Figure  34  were  increased  by  the  factor, 
practically  all  of  them  would  be  located  in  a band  above  the  curve. 

This  fact  helps  substantiate  the  validity  of  the  criterion  curve.  As 
a suggested  methodology  for  the  use  of  the  curve  in  Figure  34,  one 
could  replot  the  curve  with  values  of  H reduced  by  1.03  and  then,  using 
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the  calculated  data  of  Table  3 in  establishing  the  boundary  layer 
characteristics,  one  would  have  a curve  consistent  with  TAPS  and 
experimental  transition  information. 

Of  the  various  simplified  criteria  for  transition  developed  by 
others,  the  correlation  line  presented  by  Granville  in  Reference  15 
comes  closest  to  the  criteria  presented  here.  He  plots  a momentum 
thickness  Reynolds  number  difference  between  the  transition  and  critical 
values  against  a type  of  radius  gradient  parameter.  He  uses  no  pressure 
gradient  parameter  but  merely  collapses  all  the  data  to  a single 
curve  depending  only  on  his  radius  gradient  parameter.  It  would  appear 
that,  for  the  bodies  studied,  the  variation  in  M was  relatively  small 
for  most  of  the  transition  locations,  leading  to  the  use  of  a single 
correlation  line. 

Conclusions 

Techniques  have  been  developed  and  data  have  been  generated  that 
permit  one  to  estimate  in  a relatively  easy  manner  the  temperature 
distribution  necessary  to  stabilize  the  flow  over  an  axlsymmetric  body. 
Additionally,  these  techniques  and  data  provide  a means  for  determining 
the  local  heat  flux,  skin  friction,  the  laminar  separation  point,  and  the 
point  of  transition  for  a heated  axisymraetric  body. 

All  these  estimates  can  be  made  by  using  the  values  of  the  various 
quantities  presented  in  the  tables  and  by  following  the  interpolative 
procedures  and  criteria  that  have  been  outlined.  A simple  computer 
program  to  do  this  can  be  easily  written  or  one  can  Interpolate  between 


the  curves  that  are  also  included  in  this  memorandum. 
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TABLE  1 

Comparison  of  Quantities  Calculated  by  ARL  Code  and  Those  of  Reference  4 


Reference  4 ARL  Calculations 


•=  60.  T^  = 60) 


f"(0) 

0.33205733 

0.33205734 

n^* 

1.72078760 

1.72078700 

= 90.  T„  = 60) 

f"(0) 

0.46639971 

0.46639956 

0 "(0) 

-0.71321360 

-0.71321537 

^6* 

1.51454910 

1.51454330 

= 150.  T„  = 60) 

f"(0) 

0.74319378 

0.74319466 

e"(o) 

-0.78381642 

-0.78381826 

Tig* 

1.21124000 

1.21123790 
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TABLE  2 

Comparison  with  Blasius'  Results 

Blasius  ARL  Calculations 

5 J 1.7208  1.72078678 

fUe  1^/2 

6 0.664  0.66740880 
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X 

M 

H 

lATA 

FOR  DELTA  T = 

0.0 

1 

-0.50 

-O.OAAO 

3.55605 

2 

-0.50 

-O.OAGO 

3.2  1105 

3 

-0.5C 

-0.0300 

2.39515 

A 

-0.50 

-0.0200 

2.7  AAIA 

5 

-0.50 

-0.0100 

2.5 A908 

6 

-0.50 

0.0000 

2.58210 

7 

-0.50 

0.0500 

2. A 1156 

8 

-0.50 

0. 1000 

2.33728 

9 

-0.50 

0.2000 

2.26732 

10 

-0.50 

0.5000 

2.1 9865 

11 

-0.50 

0.7500 

2.1 7739 

12 

-C.50 

0.9000 

2.16921 

13 

0.00 

-0.0900 

3.79906 

lA 

0.00 

-0.0800 

3.20730 

15 

0.00 

-0.0700 

3.01196 

16 

0.00 

-0.0600 

2.89151 

17 

0.00 

-0.0553 

2.3 Asoa 

18 

0.00 

-0.0500 

2.80588 

19 

o.co 

-O.OAOO 

2.7A051 

20 

0.00 

-0.0300 

2.68832 

21 

0.00 

-0.0200 

2.5A538 

22 

0.00 

-0.0100 

2.60925 

23 

0.00 

0.0000 

2.57831 

2A 

0.00 

0.0500 

2.A71A8 

25 

0.00 

0.1000 

2.A0735 

26 

0.00 

0.2000 

2.33267 

27 

0.00 

0.3000 

2.28990 

28 

0.00 

0.5000 

2.2A200 

29 

0.00 

0.7500 

2.2 1033 

30 

0.00 

1.0000 

2.19130 

31 

0.50 

-0.1300 

3. A 5903 

32 

0.50 

-0.1200 

3.20AA2 

33 

0.50 

-0.1100 

3.06227 

3A 

0.50 

-0.1000 

2.96361 

35 

0.50 

-0.0900 

2.88873 

36 

0.50 

-0.0800 

2.82890 

37 

0.50 

-0.0600 

2.73771 

38 

0.50 

-0.0553 

2.72016 

39 

0.50 

-O.OAOO 

2.67031 

AO 

0.50 

-0.0200 

2.6 1779 

A1 

0.50 

0.0000 

2.575A1 

A2 

0.50 

0,0500 

2.A9752 

A3 

0.50 

0.1000 

2.AA378 

AA 

0.50 

0.2000 

2.37352 

A5 

0.50 

0.5000 

2.275A1 

q T 

(watts/ft^)  (Ib/ft^) 

0.00 

0.039180 

A. 32630 

0.00 

0.085520 

3.77832 

0.00 

0. 1537A0 

3.17978 

0.00 

0.2030A0 

2.8A531 

0.00 

0.2A3920 

2.61182 

0.00 

0.2796A0 

2.A3356 

0.00 

0.A18180 

1.90616 

0.00 

0.5232A0 

1.62523 

0.00 

0.688260 

1.30985 

0.00 

1.038010 

0.91627 

0.00 

1.257630 

0.76767 

o.co 

1.372690 

0.70708 

0.00 

0.022Ae0 

3.30322 

0.00 

0.120950 

2.67167 

0.00 

0.17A650 

2.A1957 

0.00 

0.217A20 

2.2A8AA 

0.00 

0.235320 

2.18320 

o.co 

0.25A250 

2. 1177A 

o.co 

0.287100 

2.01193 

o.co 

0.317160 

1.92315 

0.00 

0.3AA950 

1.6A683 

0.00 

0.370950 

l,78C0A 

0.00 

0.395500 

1.72079 

0.00 

0.502200 

1.A989A 

0,00 

0.591A00 

1. 3A786 

0.00 

0.7A0000 

1. 1A920 

0.00 

0.86A300 

1.01960 

0.00 

1.071500 

0.85A68 

0.00 

1.285000 

0.73010 

0.00 

1.A68000 

0.6A790 

0.00 

0.085580 

2.A2005 

0.00 

0.1A8130 

2.181A1 

0.00 

0.193990 

2.03390 

0.00 

0.232A20 

1.92A0A 

0.00 

0.266290 

1.33585 

o.co 

0.296980 

1.76201 

0.00 

0.351670 

1.6A27A 

o.co 

0.363530 

1.6187A 

0.00 

O.A00060 

1 .5A8A8 

0.00 

0.AA3920 

1.A7082 

0.00 

0.  AeA350 

1.AC502 

0.00 

0.57AA20 

1.27532 

0.00 

0.653310 

1.17777 

0.00 

0.789370 

1.0^719 

0.00 

1.10A830 

0.60A55 
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22  February  1978 
JJE;jep 


9 , T 


X 

M 

H 

(watts/ft^)  (Ib/ft^) 

46 

0.50 

0.7500 

2.2  3753 

0.00 

1.31235C 

0.69784 

47 

0.50 

0.9000 

2.22232 

0.00 

1.422630 

0.55124 

48 

1.00 

-0.0800 

2.7  3535 

0.00 

0.40608C 

1.42265 

49 

1.00 

-0.0600 

2.58313 

0.00 

0.446530 

1.35987 

50 

1.00 

-0.0553 

2.67231 

0.00 

0.458020 

1 . 34649 

51 

1.00 

-0.0400 

2.64014 

C.OO 

0.487830 

1.30591 

52 

1.00 

-0.02C0 

2.60398 

0.00 

0.52460G 

1.25868 

53 

l.OC 

0.0000 

2.57299 

0.00 

0.559300 

1.21678 

54 

1.00 

0.0500 

2.51165 

0.00 

0.638700 

1.12936 

55 

1.00 

0. 1000 

2.46566 

0.00 

0.710200 

1.05953 

56 

1.00 

0.2000 

2.40142 

0.00 

0.636400 

0.95308 

57 

1.00 

0.5000 

2.30192 

0.00 

1.137700 

0.76262 

58 

1.00 

0.7500 

2.26041 

0.00 

1.339600 

0.66971 

59 

1.00 

0.9000 

2.24327 

0.00 

1.447600 

0,62810 

60 

2.00 

-0.0800 

2.66394 

o.co 

0. 565760 

1.09494 

61 

2.00 

-0.06C0 

2.6  3615 

0.00 

0.597470 

1.06627 

62 

2.00 

-0.0400 

2.6 1137 

0.00 

0.627800 

1.04003 

63 

2.00 

-0.0200 

2.58907 

0.00 

0.656930 

1.C1587 

64 

2.00 

0.0000 

2.56890 

0.00 

0.685000 

0.99350 

65 

2.00 

0.0500 

2.5  2584 

0.00 

0.751100 

0.94401 

66 

2.00 

0.1000 

2.49077 

0.00 

0.612400 

0.90179 

67 

2.00 

0.2000 

2.43677 

0.00 

0.923900 

0.63281 

68 

2.00 

0.5000 

2.34131 

0.00 

1.20 1700 

0.69576 

69 

2.00 

0.7500 

2.29669 

0.00 

1.393400 

0.62268 

70 

2.00 

0.9000 

2.2  7734 

0.00 

1,497100 

0.56667 

71 

3.00 

-0.0800 

2.6  3278 

0.00 

0.689900 

0.92342 

72 

3.00 

-0.0600 

2.6 1393 

0.00 

0.716310 

0.90619 

73 

3.00 

-0.0400 

2.59653 

0.00 

0.741900 

0.89002 

74 

3.00 

-0.0200 

2.58043 

0.00 

0.766760 

0.87479 

75 

3.00 

0.0000 

2.56546 

0.00 

0.790900 

0.66039 

76 

3,00 

0.0500 

2.5  3224 

0.00 

0.848700 

0.62759 

77 

3.00 

0.1000 

2.50388 

0.00 

0.903200 

0.79858 

78 

3.00 

0.2000 

2.45788 

0.00 

1.004300 

0.74921 

79 

3.00 

0.5000 

2.36912 

0.00 

1.263300 

0.64A19 

60 

3.00 

0.7500 

2.32411 

0.00 

1.446000 

0.56460 

81 

3.00 

0.9000 

2.30387 

0.00 

1.545700 

0.55607 

82 

3.00 

1.1000 

2.28210 

0.00 

1.669799 

0.52390 

83 

4.00 

-0.0800 

2.6 1455 

0.00 

0.795030 

0.81353 

84 

4.00 

-0.0600 

2.60029 

0.00 

0.818130 

0.60174 

85 

4.00 

-0.0400 

2.58691 

0.00 

0.840680 

C.79C51 

66 

4.00 

-0.0200 

2.57430 

0.00 

0.662730 

0.77980 

87 

4.00 

0.0000 

2.56241 

0.00 

0.66430C 

0.76956 

68 

4.00 

0.0500 

2.5  3539 

0.00 

0.935300 

0. 7‘t579 

89 

4.00 

0.1000 

2.5 1164 

0.00 

0.985900 

C. 72429 

90 

4.00 

0.2000 

2.47171 

0.00 

1.079000 

0.66670 

91 

4.00 

0.5000 

2.38969 

0.00 

1.3220CC 

0.60276 

92 

4.00 

0.7500 

2.34551 

o.co 

1 .497000 

0. 55292 

93 

4.00 

0.9000 

2.32507 

0.00 

1.593000 

0.52849 
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22  February  1978 
JJE:jep 


X 

M 

H 

‘I  2 

(watts/ft  ) 

T 

1 (Ib/ft^) 

rig* 

94 

4.00 

0.9600 

2.3 1782 

0.00 

1.630420 

0.51961 

95 

4.00 

1.1200 

2.30062 

0.00 

1.725300 

0.49600 

96 

4.00 

1.2800 

2.28592 

0.00 

1.815420 

0.47891 

97 

5.00 

-0.0800 

2.60227 

0.00 

0.867850 

0.73541 

98 

5.00 

-0.0600 

2.59082 

0.00 

0.908640 

0.72670 

99 

5.00 

-0.0400 

2.57993 

0.00 

0.929030 

0.71633 

100 

5.00 

-0.0200 

2.56957 

0.00 

0.949050 

0.71027 

101 

5.00 

0.0000 

2.55970 

0.00 

0.968700 

0.70251 

102 

5.00 

0.0500 

2.5  3691 

0.00 

1.016400 

0.68428 

103 

5.00 

0.1000 

2.5 1646 

0.00 

1.062100 

0.66752 

104 

5.00 

0.2000 

2.46125 

0.00 

1.146800 

0.63766 

105 

5.00 

0.5000 

2.40543 

0.00 

1.379300 

0.56854 

106 

5.00 

0.7500 

2.36263 

0.00 

1.547200 

0.52598 

107 

5.00 

0.9000 

2.34235 

0.00 

1.640200 

0.50474 

108 

5.00 

0.9600 

2.33510 

0.00 

1.675980 

0.49696 

109 

5.00 

1.12UO 

2.31773 

0.00 

1.768150 

0.47789 

110 

5.00 

1.2800 

2.30270 

0.00 

1.855950 

0.46090 

111 

5.00 

1.4400 

2.28958 

0.00 

1.939940 

0 . 44  564 

DATA  FOR  DELTA  T = 


112 

-0.50 

-0.0400 

113 

-0.50 

-0.0200 

114 

-0.50 

0.0000 

115 

-0.50 

0.0500 

116 

-0.50 

0.1000 

117 

-0.50 

0.2000 

118 

-0.50 

0.5000 

119 

0.00 

-0.0900 

120 

0.00 

-0.0800 

121 

0.00 

-0.0700 

122 

0.00 

-0.0600 

123 

0.00 

-0.0553 

124 

0.00 

-0.0400 

125 

0.00 

-0.0200 

126 

0.00 

0.0000 

127 

0.00 

0.0500 

128 

0.00 

0.1000 

129 

0.00 

0.2000 

130 

0.00 

0.5000 

131 

0.50 

-0.  1200 

132 

0.50 

-0.  1000 

133 

0.50 

-0.C900 

134 

0.50 

-0.08GO 

135 

0.50 

-0.0600 

136 

0.50 

-0.0553 

137 

0.50 

-0.0400 

138 

0.50 

-0.0200 

139 

0.50 

0.0000 

0.0 

2.37276  1237.00 

2.52498  1469.71 

2.39277  1605.00 

2.25160  1834.00 

2.19045  2003.99 

2.13379  2273.00 

2.07987  2876.  19 

3.13521  1570.00 

2.86927  1749.00 

2.73210  1860.00 

2.64085  1946.14 

2.60701  1980.93 

2.52147  2078.50 

2.44A21  2182.52 

2.38908  2270.00 

2.30061  2449.00 

2.24743  2594.00 

2.18581  2834.00 

2.11235  3379.00 

2.86660  2142.00 

2.69544  2316.00 

2.63620  2364.00 

2.59147  2443.00 

2.5  1677  2545.63 

2.50461  2567.35 

2.46414  2633.27 

2.42115  2710.62 

2.38626  2780.00 


0.100700 

3.25866 

0.199140 

2.50159 

0.265550 

2.14189 

0.386130 

1.67379 

0.477550 

1.42349 

0.621040 

1.14311 

0.925300 

0.79553 

0.081500 

2.63566 

0.142420 

2.30422 

0.185880 

2.11028 

0.221870 

1.96995 

0.237110 

1.91517 

0.281630 

1.76888 

0.331660 

1.62532 

0.375500 

1.51455 

0.468400 

1.31755 

0.546100 

1.16355 

0.675300 

1.00655 

0.963650 

0.74492 

0. 174420 

1.88139 

0.243130 

1.66129 

0.271730 

1.60646 

0.297900 

1.54633 

0.344920 

1.44429 

0.355160 

1.42355 

0.366760 

1.36250 

0.424810 

1.29454 

0.459940 

1.23662 
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22  February  1978 
JJE: jep 


X 

M 

H 

’ 2 

(watts/ft  ) 

(Ib/ft^) 

Hg* 

1^0 

0.50 

0.0500 

2.32177 

2932.00 

0.538310 

1 . 12185 

lAl 

0.50 

0. 1000 

2.27714 

3062.00 

0.606990 

1.03515 

1A2 

0.50 

0.2000 

2.2 1888 

3282.00 

0.72542C 

0.90993 

1^3 

0.50 

0.5000 

2.1 3835 

3793. 19 

0.999780 

0.70265 

14A 

1.00 

-0.0800 

2.5 1652 

2939.39 

0. 398280 

1.25079 

145 

1.00 

-0.06C0 

2.47423 

3016.60 

0.43497C 

1.19638 

146 

l.OC 

-0.0553 

2.46542 

3033.64 

0.443190 

1 .13473 

147 

1.00 

-0.0400 

2.43915 

3086.60 

0.469040 

1.14928 

146 

1.00 

-0.02C0 

2.409-^3 

3150.79 

0.500960 

1.10784 

149 

1.00 

0.0000 

2.38387 

3211.00 

0.531100 

1.07095 

150 

1.00 

0.05C0 

2.33310 

3345.00 

0.600200 

0.99363 

151 

1.00 

0.1000 

2.29508 

3463.00 

0.662400 

0.93165 

152 

1.00 

0.20C0 

2.24154 

3668.00 

0.772300 

0.63689 

153 

1.00 

0.5000 

2.1 5933 

4154.70 

1.034420 

0.66724 

154 

l.CO 

0.9000 

2.1 1165 

4664.75 

1.3C386C 

0.54778 

155 

2.00 

-0.0800 

2.45799 

3724.00 

0.546960 

0.96343 

156 

2.00 

-0.0600 

2.43530 

3780.23 

0.57445C 

0.95838 

157 

2.00 

-0.0400 

2.41491 

3833.39 

0.600780 

0.91538 

158 

2.00 

-0.0200 

2.39656 

3883.89 

0.62608C 

0.89413 

159 

2.00 

0.0000 

2.37991 

3932.00 

0.650500 

0.87442 

160 

2.00 

0.0500 

2.34424 

4044.00 

0.70790C 

0.83070 

161 

2.00 

0.1000 

2.3 1512 

4146.00 

0.761300 

0.79326 

162 

2.00 

0.2000 

2.2  7015 

4330.00 

0.858400 

0.73196 

163 

2.00 

0.5000 

2.19093 

4778.79 

1.100160 

0.60987 

164 

3.00 

0.0000 

2.37650 

4540.00 

0.751100 

0.75727 

165 

3.00 

0.0500 

2.34901 

4639.00 

0.601400 

0.72830 

166 

3.00 

0.1000 

2.32548 

4730.00 

0.848800 

0.70260 

167 

3.00 

0.2000 

2.28722 

4697.00 

0.936800 

0.65877 

168 

3.00 

0.5000 

2.2 1333 

5317.79 

1.16218C 

0.56528 

169 

4.00 

0.0000 

2.37358 

5076.00 

0.83970G 

0.67733 

170 

4.00 

0.0500 

2.35119 

5165.00 

0.885COO 

0.65634 

171 

4.00 

0.1000 

2.3  3148 

5248.00 

0.928100 

0.63731 

172 

4.00 

0.2000 

2.29826 

5403.00 

1 .C0900C 

0.60396 

173 

4.00 

0.5000 

2.22991 

5800.29 

1.221060 

0.52930 

174 

5.00 

0.0000 

2.37093 

5561.00 

0.919900 

0.61831 

175 

5.00 

0.0500 

2.35206 

5642.00 

0.96130C 

0.60221 

176 

5.00 

0.1000 

2.33509 

5719.00 

1.001200 

0.58739 

177 

5.00 

0.20C0 

2.30578 

5664.00 

1.0766CC 

0.56092 

178 

5.00 

0.5000 

2.24256 

6241.59 

1.277270 

0.49946 

179 

5.00 

0.9000 

2.1 8999 

6662.60 

1.504290 

0.44262 

DATA 

FOR  DE 

LTA  T = 

60.0 

160 

-0.50 

-0.0400 

2.S  3442 

2756.00 

0. 108640 

2.85912 

181 

-0.50 

-0.0200 

2.36421 

3200.68 

0.192570 

2.22318 

182 

-0.50 

0 . COCO 

2.25545 

3473.00 

0.250550 

1 . 90540 

183 

-0.50 

0.0500 

2.1 3895 

3941.00 

0.356270 

1.48659 
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JJE: jep 

q T 


X 

M 

H 

2 

(watts/ft  ) 

(Ib/ft^) 

rig* 

164 

-0.50 

0.1000 

2.08935 

4290.00 

0.436470 

1.26198 

185 

-0.50 

0.2000 

2.04470 

4851.00 

0.562420 

1.01080 

186 

-0.50 

0.5000 

2.00433 

6117.92 

0.829970 

0.70095 

187 

0.00 

-0.09C0 

2.30413 

3603.00 

0.106700 

2.26039 

188 

0.00 

-0.0800 

2.63106 

3897.00 

0.153640 

2.02170 

189 

0.00 

-0.0700 

2.52847 

4103.00 

0.189900 

1.86443 

190 

0.00 

-0.0600 

2.45699 

4266.68 

0.220630 

1 . 74603 

191 

0.00 

-0.0553 

2.42999 

4334.33 

0.233750 

1.69906 

192 

0.00 

-0.0400 

2.36075 

4526.45 

0.272330 

1.57202 

193 

0.00 

-0.0200 

2.29735 

4734.30 

0.315960 

1.44564 

194 

0.00 

0.0000 

2.25178 

4911.00 

0.354300 

1.34731 

195 

0.00 

0.0500 

2.1 7853 

5275.00 

0.435700 

1.17133 

196 

0.00 

0.1000 

2.1 3470 

5573.00 

0.503800 

1.05118 

197 

0.00 

0.2000 

2.08457 

6067.00 

0.617300 

0.89235 

198 

0.00 

0.5000 

2.02663 

7203.29 

0.870370 

0.65811 

199 

0.50 

-0. 12C0 

2.62851 

4773.00 

0.168170 

1.65071 

200 

0.50 

-0.1000 

2.49952 

5096.00 

0.245720 

1.48738 

201 

0.50 

-0.0900 

2.45440 

5226.00 

0.270210 

1.42563 

202 

0.50 

-0.0800 

2.4 1703 

5341.00 

0.292800 

1.37228 

203 

0.50 

-0.0600 

2.35809 

5543.75 

0.333540 

1.28355 

204 

0.50 

-0.0553 

2.34653 

5586.94 

0.342440 

1.26538 

205 

0.50 

-0.0400 

2.3 1330 

5718.54 

0.369990 

1.21166 

206 

0.50 

-0.0200 

2.27785 

5873.96 

0.403240 

1.15154 

207 

0.50 

0.0000 

2.24696 

6015.00 

0.433970 

1.10008 

208 

0.50 

0.0500 

2.19549 

6323.00 

0.502630 

0.99763 

209 

0.50 

0.1000 

2.1 5857 

6586.00 

0.562850 

0.91998 

210 

0.50 

0.2000 

2.1 1068 

7042.00 

0.666740 

0.80765 

211 

0.50 

0.5000 

2.04579 

8102.00 

0.907530 

0.62195 

212 

1.00 

-0.0800 

2.35586 

6401.40 

0.385140 

1.11159 

213 

1.00 

-0.0600 

2.32121 

6555.09 

0.417090 

1.06381 

214 

1.00 

-0.0553 

2.3 1396 

6589.18 

0.424270 

1.05354 

215 

1.00 

-0.0400 

2.29231 

6695.29 

0.446830 

1.02222 

216 

1.00 

-0.0200 

2.26775 

6824.79 

0.474740 

0.98550 

217 

1.00 

0.0000 

2.24659 

6945.00 

0.501100 

0.95269 

218 

1.00 

0.C500 

2.20447 

7218.00 

0.561600 

0.86372 

219 

1.00 

0.1000 

2.1 7293 

7459.00 

0.616200 

0.82826 

220 

l.CO 

0.2000 

2.12871 

7881.00 

0.712500 

0.74330 

221 

1.00 

0.5000 

2.0  6178 

8886.96 

0.942550 

0.59108 

222 

1.00 

0.9000 

2.02397 

9949.09 

1.179110 

0.46414 

223 

2.00 

-0.08C0 

2.30720 

8087.29 

0.523250 

0.85678 

224 

2.00 

-0.060C 

2.28845 

8200.10 

0.547250 

0.83464 

225 

2.00 

-0.0400 

2.27163 

8307.09 

0.570260 

0.81426 

226 

2.00 

-0.0200 

2.2  5643 

8409.00 

0.592390 

0.79540 

227 

2.00 

0.0000 

2.24261 

6506.00 

0.613700 

0.77787 

228 

2.00 

0.0500 

2.2 1299 

8734.00 

0.664100 

0.73888 

229 

2.00 

0.1000 

2.1 8876 

8942.00 

0.710800 

0.70543 

230 

2.0C 

0.2000 

2.1 5152 

9317.00 

0.796000 

0.65052 

231 

2.00 

0.5000 

2.08629 

10242.50 

1.006090 

0.54099 
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X M H 


232 

3.00 

0.0000 

2.2  3922 

233 

3.00 

0.0500 

2.2 1633 

23A 

3.00 

0. 1000 

2.1 9683 

235 

3.00 

0.2000 

2.16503 

236 

3.00 

0.5000 

2.10388 

237 

4.00 

0,0000 

2.2  3631 

238 

4.00 

0.05C0 

2.2 1770 

239 

4.00 

0.1000 

2.20130 

2^0 

4.00 

0.20C0 

2.17366 

241 

4.00 

0.5000 

2.1 1693 

242 

5.00 

0.0000 

2.23367 

243 

5.00 

0.0500 

2.2 1797 

244 

5.00 

0, 1000 

2.20365 

245 

5.00 

0.2000 

2.1 7945 

246 

5.00 

0.5000 

2.1 2686 

247 

5.00 

0.9000 

2.08338 

DATA 

FOR  DELTA  T = 

90.0 

248 

-0.50 

-0.0400 

2.46339 

249 

-0. 5C 

-0.0200 

2.24701 

250 

-0.50 

0.0000 

2.1 5691 

251 

-0.50 

0.0500 

2.06114 

252 

-0.50 

0.1000 

2.02165 

253 

-0.50 

0.2000 

1.98775 

254 

0.00 

-0.0900 

2.58584 

255 

0.00 

-0.0800 

2.46010 

256 

0.00 

-0.C700 

2.37993 

257 

0.00 

-0.0600 

2.32242 

258 

0.00 

-0.0553 

2.30043 

259 

0.00 

-0.04C0 

2.24353 

260 

0.00 

-0.0200 

2.19101 

261 

0.00 

0 . OCuO 

2.1 5317 

262 

o.co 

0.0500 

2.09258 

263 

0.00 

0.1000 

2.0  5676 

264 

0.00 

0.2000 

2.0 1666 

265 

0.00 

0.5000 

1.97261 

266 

0,50 

-0. 12C0 

2.45756 

267 

0.50 

-0.1000 

2.35628 

268 

0.50 

-0.0900 

2.3 1985 

269 

0.50 

-0.0800 

2.28935 

270 

0.50 

-0.0600 

2.24086 

271 

0.50 

-0.0553 

2.23130 

272 

0.50 

-0.0400 

2.20375 

273 

0.50 

-0.0200 

2.1 7431 

274 

0.50 

O.OOCO 

2.15030 

275 

0.50 

C.C500 

2.10597 

9 . T 


2 

(watts/ft  ) 

(Ib/ft^) 

9822.00 

0.708700 

0.67366 

10022.00 

0.752700 

0.64783 

10206.00 

0.79430C 

0.62438 

10549.00 

0.671400 

0.56567 

11413.30 

1.06911C 

0.50181 

10982. CO 

0.79230C 

0.60254 

11161.00 

0.83190C 

0.58384 

11331.00 

0.669700 

0.56665 

11646.00 

0.940700 

0.53703 

12461.19 

1.126610 

0.47010 

12030.00 

0.867900 

0.55004 

12194.00 

0.904200 

0.53570 

12351.00 

0.939100 

0.52247 

12646. CO 

1 .C053C: 

0.49881 

13420.00 

1.181210 

0.44375 

14287.59 

1.380370 

0.39274 

4500.00 

0.111780 

2.54460 

5141.33 

0. 18441C 

1.99676 

5546.00 

0.235480 

1.71296 

6259.00 

0.32698C 

1.33522 

6796.00 

0.40004C 

1.13212 

7661.00 

0.511790 

0.90528 

5976.00 

0.119670 

1.98775 

6364.00 

0.156080 

1.79930 

6653.00 

0.189000 

1.66708 

6889.29 

0.215590 

1.56490 

6987.94 

0.227030 

1.52388 

7270.94 

0.260800 

1.41192 

7580.66 

0.29917C 

1.29937 

7846. CO 

0.333000 

1.21124 

8397.00 

0.404900 

1.05271 

8851.00 

0.465200 

0.94414 

9610.00 

0.565700 

0.80053 

11365.39 

0.790390 

C. 58906 

7794. CO 

0. 19361C 

1.46913 

8250.00 

0.242810 

1.33120 

8438.00 

0.264000 

1.27773 

8607.00 

0.283700 

1.23110 

8905.03 

0.31941C 

1.15283 

8969.09 

0.327230 

1.13671 

9165.09 

0.351460 

1.08690 

9397.77 

0.380750 

1.03514 

9610. CO 

C. 40780 C 

0.98897 

10076.00 

0.468530 

C. 89675 
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w T 


X 

M 

H 

(watts/ft  ) 

(Ib/ft^) 

^6* 

276 

0.50 

0.1000 

2.07555 

10479.00 

0.521810 

0.82665 

277 

0.50 

0.2000 

2.3  3658 

11172.00 

0.613800 

0.72510 

278 

0.50 

0.50C0 

1.98556 

12604.60 

0.627310 

0.55728 

2 79 

1.00 

-0.0600 

2.23861 

10282.70 

0. 36882C 

0.99838 

260 

1.00 

-0.06C0 

2.2  0991 

10511.20 

0.396910 

0.95592 

281 

1.00 

-0.0553 

2.20369 

10561.99 

0.403220 

0.94676 

282 

1.00 

-0.0400 

2.18589 

10720.70 

0.423090 

0.91860 

283 

1.00 

-0.0200 

2.16549 

10914.79 

0.447690 

0.88591 

284 

1.00 

0.0000 

2.14789 

11096.00 

0.471000 

0.85647 

285 

1.00 

0.0500 

2.1 1291 

11508.00 

0.524400 

0.79441 

286 

1.00 

0.1000 

2.08683 

11875.00 

0.572700 

0.74438 

267 

1.00 

0.2000 

2.05053 

12518.00 

0.658000 

0.66761 

288 

1.00 

0.5000 

1.99699 

14062.50 

0.861810 

0.52999 

289 

2.00 

-0.0600 

2.18199 

13130.09 

0.518170 

0.75019 

290 

2.00 

-0.0400 

2.1 6799 

13290.40 

0.538460 

0.73195 

291 

2.00 

-0.0200 

2.1 5532 

13443.60 

0.557970 

0.71505 

292 

2.00 

0.0000 

2.14381 

13590.00 

0.576600 

0.69931 

293 

2.00 

0.0500 

2.1 1916 

13934.00 

0.621300 

0.66424 

294 

2.00 

0.1000 

2.09908 

14249.00 

0.662600 

0.63410 

295 

2.00 

0.2000 

2.06830 

14820.00 

0.738000 

0.58453 

296 

2.00 

0.5000 

2.01516 

16235.39 

0.925790 

0.48550 

297 

3.00 

0.0000 

2.14041 

15693.00 

0.666000 

0.60562 

298 

3.00 

0.0500 

2.12138 

15993.00 

0.704900 

0.56240 

299 

3.00 

0.1000 

2.1 0600 

16267.00 

0.740600 

0.56240 

300 

3.00 

0.2000 

2.0  7873 

16793.00 

0.809900 

0.52635 

301 

3.00 

0.5000 

2.02847 

18112.39 

0.98491C 

0.45057 

302 

4.00 

0.0000 

2.13744 

17545.00 

0.744700 

0.54168 

303 

4.00 

0.0500 

2.12192 

17816.00 

0.779600 

0.52487 

304 

4.00 

0.1000 

2.10626 

18072.00 

0.813000 

0.50958 

305 

4.00 

0.2000 

2.08528 

18551.00 

0.875800 

0.48270 

306 

4.00 

0.5000 

2.03839 

19792.39 

1.040320 

0.42223 

307 

5.00 

0.0000 

2.1 3477 

19220.00 

0.815700 

0.49449 

308 

5.00 

0.0500 

2.1 2166 

19468.00 

0.647800 

0.46159 

309 

5.00 

0.1000 

2.10988 

19705.00 

0.878600 

0.46969 

310 

5.00 

0.2000 

2.08955 

20152.00 

0.937100 

0.44637 

311 

5.00 

0.5000 

2.34595 

21328.79 

1.092750 

0.39865 
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TABLE  4 

Laminar  Boundary  Layer  Data  for  Various  Ambient  Temperatures 


AT 

M 

H 

(watts/ft  ) 

(Ib/ft^) 

rig* 

DATA 

FCR  LAFBCA  = 0. 

0,  T AT 

INFINITY 

= 60.0 

1 

C.O 

0.0000 

2.57631 

O.CC 

0.395500 

1.72079 

2 

0.0 

0.0500 

2.47148 

C.CO 

0.5C22C0 

1.49841 

3 

0.0 

O.IOCO 

2.40725 

C.CO 

0.591400 

1.34786 

A 

0.0 

0.2000 

2.33267 

0.00 

0.740000 

1.14920 

5 

30.0 

0.0000 

2.389C8 

2270.00 

0.375500 

1.51455 

6 

30.0 

0.0500 

2.30061 

2449.00 

0.468400 

1.31755 

7 

30.0 

O.ICOO 

2.24743 

2594.00 

0.546100 

1.18355 

6 

30.0 

0.2000 

2.1 8561 

2834.00 

0.675300 

1.00655 

9 

60.0 

c.ocoo 

2.2  5178 

4911.00 

0.354300 

1.34731 

1C 

60.0 

0.0500 

2.1 7852 

5275. CO 

0.435700 

1.17133 

11 

60.0 

O.IOCO 

2.13470 

5573.00 

0.503800 

1.05118 

12 

60.0 

0.2000 

2.D  8457 

6067.00 

0.617300 

0.89235 

13 

90.0 

c.oooo 

2.15317 

7646.00 

C.3330C0 

1.21124 

14 

90.0 

0.0500 

2.09258 

8397.00 

0.404900 

1.05271 

15 

90.0 

0.1000 

2.0  5676 

8851.00 

0.465200 

0.94414 

16 

90. C 

C.20C0 

2.01666 

9610.00 

0.565700 

0.80053 

DATA 

FOR  LAFBCA  = 0. 

0,  T AT 

INFINITY 

= 90.0 

17 

0.0 

0.0000 

2.57821 

O.CC 

0.325200 

1.72079 

18 

C.O 

0.0500 

2.47148 

0.00 

0.412900 

1.49841 

19 

0.0 

0.1000 

2.40735 

0.00 

0.486300 

1.34786 

20 

C.O 

0.2000 

2.33267 

0.00 

0.608400 

1.14920 

21 

30.0 

0.0000 

2.43051 

2429.00 

0.309600 

1.54479 

22 

30.0 

0.0500 

2.33984 

2617.00 

0. 387500 

1.34442 

23 

30.0 

0.1000 

2.26544 

2769.00 

0.452700 

1.20821 

24 

30.0 

0.2000 

2.2  225  1 

3022.00 

0.561300 

1.02828 

25 

60.0 

C.OOOO 

2.32214 

5165.00 

0.293600 

1.40CC9 

26 

60.0 

0.0500 

2.24502 

5544.00 

0.363300 

1.21811 

27 

60.0 

0. 1000 

2.19901 

5855.00 

0.421700 

1.09401 

28 

60.0 

0.2000 

2.14645 

6371.00 

0.519000 

0.92994 

29 

90.0 

C.OOOO 

2.24388 

8141.99 

0.277800 

1.28059 

3C 

90.0 

0.0500 

2.17820 

8714.00 

0.340700 

1.11406 

31 

90.0 

O.IOCO 

2.13947 

9187.00 

0.393400 

1 .00018 

32 

90.0 

0.2000 

2.09609 

9974.00 

0.481300 

0.84953 

DATA 

FCR  LAKBCA  = 0. 

Ot  T AT 

INFINITY 

= 120.0 

33 

C.C 

C.OCCG 

2.57831 

0.00 

0.2773CC 

1 .72079 

34 

0.0 

0.05CC 

2.47140 

O.CC 

C.3520CC 

1.49841 

35 

C.O 

0.1000 

2.40735 

0.00 

0.414600 

1.34786 
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AT 

M 

H 

36 

C.O 

C.20C0 

2.33267 

37 

3C,0 

c.cooo 

2.46277 

38 

30.0 

C.05C0 

2.37008 

39 

30.0 

0.1000 

2.3 1457 

40 

30.0 

C.20CC 

2.25C4C 

41 

60.0 

o.ooco 

2.37774 

42 

60.0 

c.icoo 

2.24912 

43 

60. 0 

0.2000 

2.1 9435 

44 

90. 0 

c.oooo 

2.3 1655 

45 

90.0 

0. 1000 

2.20496 

46 

90.0 

C.20C0 

2.1 5871 

2 ^ 
^ t rM.^\  /ii_/ 


2 

(watts/ft  ) 

(Ib/ft^) 

C.CC 

0.518800 

1.1‘.920 

2565.00 

0.264700 

1.56982 

2759.00 

0.332300 

1.36651 

2918.00 

0.388800 

1.22840 

3180.00 

0.482900 

1.04600 

5383.00 

0.25220C 

1.44408 

6096.00 

0.3650CC 

1.12939 

6629.00 

0.450700 

0.96086 

8400.00 

0.239900 

1.33893 

9475.00 

0.343300 

1.04695 

10284.00 

0.421900 

0.89027 

i 


-35- 


22  February  1978 
JJE: jep 


TABLE  5 

Critical  Reynolds  Number  Variation  with  Shape  Factor 


1 

105.000 

3.80000 

2 

115.000 

3.50000 

3 

130.000 

3.20000 

150.000 

3.00000 

5 

180.000 

2.90000 

6 

232.000 

2.80000 

7 

326.000 

2.70000 

6 

A85.000 

2.6000  0 

9 

675.000 

2.55000 

10 

1000.000 

2.51800 

11 

1500,000 

2.48750 

12 

2000,000 

2.46600 

13 

3000.000 

2.43750 

lA 

AOOO.OOO 

2.41800 

15 

5000.000 

2.40250 

16 

6000.000 

2.3B950 

17 

7000.000 

2.37550 

le 

8000.000 

2.36050 

19 

9000.000 

2.34400 

20 

10000.000 

2.32540 

21 

10500,000 

2.31500 

22 

11000.000 

2.29750 

23 

11210.000 

2.26000 
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TABLE  6 

Transition  Reynolds  Number  Variation  with  Shape  Factor 

H Rx.  logTn[R„.  ] 

*trans  “10  *trans 


2.20CC0 

e. 51732 

0.32909E 

09 

2.25C00 

e.2eoA5 

0. 1907AE 

09 

2.30CCC 

8.0AA05 

C.11C67E 

09 

2.35CCC 

7.80882 

0.64390E 

08 

2.ACCCC 

7.575A6 

0.37625E 

08 

2.A5CC0 

7.3AA73 

C.22117E 

08 

2.50C0C 

7.11730 

0.131CCE 

08 

2.55CC0 

6.69390 

0.78324E 

07 

2.600C0 

6.6752A 

0.47341E 

07 

2.65CCC 

6.A620A 

0.28976E 

07 

2.70CC0 

6.25501 

C.  17989E 

07 

2.75000 

6.05486 

0. 11346E 

07 

2.8CCC0 

5.66231 

0.72e29E 

06 

2.eA999 

5.67808 

C.47651E 

C6 

2.90CCC 

5.50287 

0.318326 

06 

2.95C00 

5.33741 

0.21747E 

06 

Figure  .1  - Shape  Factor  Variation  with  M for  AT=0 


Factor  Variation  with  M for  AT=30‘’F 


Factor  Variation  wit 


Factor  \Jarlation  with  M for  AT=90°F 


Friction  Variation  with  M for  AT“0 


Friction  Variation  with  M for  AT“30®F 


u.o 


ebruary  1978 
jep 

3> 


Local  Skin  Friction  Variation  with  M for 


Figure  11  - Local  Skin  Friction  Variation  with  M for  AT=90®F 


Variation  with  M for  AT=60°F 


Factor  Variation  with  M for  T =60°F 


FALKNER-SKAN  PARAMETER, 


Factor  Variation  with  M for  T “120“F 


50.0  fl/sec 


Flux  Variation  with  M for  T "60®F 
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ure  22  - Local  Skin  Friction  Variation  with  M for  T *60*F 


Friction  Variation  with  M for  T -90“? 


Friction  Variation  with  M for  T “120‘’F 
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Figure  28  - Correlation  of  Shape  Factor  with  Critical  Reynolds  Number 


Interpolation  Procedure 


Figure  33  - Transition  Reynolds  Number  Variation  with  H from  e Stability 
Calculations 


BODY  NO.  1 
BODY  NO.  2 

MODIFIED  SEARS-HAACK  BODY  ] 
MODIFIED  PARABOLIC  BODY  > REF.  12 
MODIFIED  ELLIPTIC  BODY 


Figure  34  - Log  Plot  of  Transition  Reynolds  Number  Variation  with  Shape  Factor 


60'^F 


Transition  Reynolds  Number  Information 


Transition  Reynolds  Number  Information 


Transition  Reynolds  Number  Information  for  AT»60“F 


Transition  Reynolds  Number  Information  for  AT«60®F 


RADIUS  GRADIENT  PARAMETER, 


i 


1 

I 

I 


r 
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APPENDIX  A 

Property  Variation  with  Temperture  for  Water 


Without  getting  into  all  the  details  of  why,  the  formulation  of 
the  variation  of  the  water  properties  with  temperature  by  Touloukian 
and  others  (References  9 and  10)  was  used  in  Reference  4.  It  was 
judged  the  "best"  from  the  standpoint  of  derivatives  of  these  properties. 
Actually,  for  the  free-stream  properties,  the  f emulations  by  Kaups  and 
Smith  (Reference  11)  were  used.  The  option  to  use  either  of  the  above 
within  the  boundary  layer  is  open  in  the  program,  but  the  equations  of 
Touloukian  and  Makita  were  the  ones  used.  The  following  are  these 
formulations: 

Specific  Heat, 


C = 2.13974  - 9.68137  x 10“^  T + 2.68536  x lo  ^ T^ 

P 

- 2.42139  X 10"®  T®  v-<--A) 


C is  in  cal.  g ^ ®K 
P ® 

T is  in  ®K 


Density,  p 


p “ 1 


(T-3.9863)^  (T+288.9414) 
508929.2  (T+68. 12963) 


0.011445 


exp 


374.3 


(2-A) 


p is  in  g/ml 


T is  in  'C 


Dynamic  Viscosity,  y 


log  VI  = - 1.64779  + 


Ihl.yi 

T-133.98 


(3-A) 


y Is  In  cp 


T Is  In  "K 


Thermal  Conductivity,  k 


k = - 9.901090  + 0.1001982  T - 1.873892  x lo"^ 


+ 1.039570  X lo"^ 


(4-A) 


k Is  In  m watts  cm  ^ “K  ^ 


T Is  In  “K 


Unit  Conversions 

The  following  conversions  are  made  so  that  the  end  quantities  are 
in  commonly  used  units: 


p X 1.9409  to  get  p in  slugs/ft^ 
y X (2.088  X lO  to  get  y in  Ib/sec/ft^ 


k X (0.16933)  to  get  k In  watts/ft/^R 


1 


1 


I 


i 
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